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Abstract. We study how gravitational focusing (GF) of dark matter by the Sun affects the 
annual and biannual modulation of the expected signal in non-directional direct dark matter 
searches, in the presence of dark matter substructure in the local dark halo. We consider the 
Sagittarius stream and a possible dark disk, and show that GF suppresses some, but not all, 
of the distinguishing features that would characterize substructure of the dark halo were GF 
neglected. 
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1 Introduction 


Astrophysical and cosmological observations indicate that dark matter (DM) is the dominant 
form of matter in the Universe. One of the most well-motivated candidates for DM is the 
weakly interacting massive particle (WIMP) [T]. WIMPs are particles with weak-scale inter¬ 
action cross sections, and with masses roughly between a few GeV and hundreds of TeV. They 
could be detected through their scattering with atomic nuclei using sensitive, low threshold 
detectors. Many such direct detection experiments are currently in operation, employing a 
variety of target nuclei and detection techniques, attempting to gain further insight into the 
exact nature of DM. 

Due to Earth’s rotation around the Sun, we expect the DM flux seen at Earth, and 
therefore the scattering rate at non-directional direct detection experiments, to be annually 
modulated. While the signal of WIMP scattering off of target nuclei can be faked, for example 
by scattering of neutrons emitted by radioactive processes in the vicinity of the detector, an 
annual modulation in the rate with the expected features is a much more difficult signature 
to be reproduced by spurious sources. Moreover, while the energy spectrum of DM events 
depends on the WIMP mass and interactions, a modulation in the rate will be present 
regardless of these details and it is therefore sometimes claimed to be a model-independent 
signature of DM. 

In order to determine the modulation amplitude and its spectrum in energy measured 
by a particular experiment, a model for the DM halo must be assumed. The standard choice 
for the main virialized component is the Standard Halo Model (SHM), an isothermal sphere 
at rest with respect to the galaxy with an isotropic Maxwellian velocity distribution [2]. 
Despite its simple form, the SHM is believed to capture the relevant characteristics of the 
dark halo and thus has been widely used in the literature. However, we expect the actual halo 
to deviate from this simple model. The local density and velocity distribution could actually 
be very different if Earth is within a DM clump (although this is unlikely [3]), stream, dark 
disk (DD), and/or tidal debris [3HT8]. 

Both a stream and a DD are well motivated candidates for DM velocity substructure, 
capable of significantly altering the expected modulation [3Itl4]. Observations of the Sagit¬ 
tarius (Sgr) dwarf galaxy show that tidally stripped stars are currently passing through the 
galactic plane not far from the Sun. Simulations of this merger suggest that the DM com¬ 
ponent of the Sgr stream may be passing through the Solar System and could contribute as 
much as 5% to the local DM density nzms]. A DD is a subcomponent of the halo that 
has a spatial distribution roughly coincident with the visible disk, co-rotating with it, but 
with a lagging angular velocity dnHm. Numerical simulations have shown that DDs can 
form in Milky Way-type galaxies from mergers of satellite galaxies |llH13j (although recent 
measurements suggest this may be unlikely |20l I21| ). Alternatively, if a subdominant portion 
of DM is dissipative in nature it has the potential to collapse and form a DD, a process 
comparable to the formation of the baryonic disk [14| . 

Ref. [22] recently performed a Eourier analysis of the expected modulation and con¬ 
sidered how the annual and higher harmonics are influenced by the eccentricity of Earth’s 
orbit and the possible existence of DM velocity substructure. Since the eccentricity is small, 
e ~ 0.016722, it is not expected to impact the leading harmonic. DM substructure, however, 
was shown in Refs. IZIESI to profoundly impact all harmonics, including the unmodulated 
rate. Furthermore, Ref. |22j pointed out that if the DM velocity distribution is smooth and 
isotropic in the galactic frame, there exist ratios of the amplitudes of the harmonics that 


- 2 - 


are independent of the scattering energy, and thus they concluded that these ratios could be 
used to probe the level and nature of anisotropy in the DM halo. 

The annual modulation of a DM signal can also be affected by the gravitational focusing 
(GF) of DM by the Sun |23H25j . The extent to which the annual (first) and biannual (second) 
harmonics in the SHM are influenced by GF was studied in Refs. |24| I25| . They found that 
GF has a nearly negligible effect on the amplitude of the first harmonic, but can signihcantly 
enhance the amplitude of the second harmonic and generate energy-dependent phases in both 
the hrst and second harmonics, especially at low energy. 

While the Fourier analysis of the rate has been studied for the SHM by taking into 
account both GF and the eccentricity of Earth’s orbit |24[I25| . the modihcation of the leading 
harmonics in the presence of DM substructure has only been studied for an eccentric orbit 
without GF [22]. The aim of this paper is to study the effect of GF on a DM halo hosting 
substructure. We begin by considering the relative importance of GF and the eccentricity of 
Earth’s orbit for the first two harmonics within the SHM, and analyze the extent to which 
the conclusions of Ref. [22| hold when GF is considered. We then study the effect of the Sgr 
stream and a DD on the annual and biannual harmonics. 

In Section we review the Fourier expansion of the rate and the procedure by which 
the effect of GF is incorporated. Section presents the amplitudes and phases of the first 
and second harmonics for the SHM with and without DM substructure. We specifically 
discuss how GF affects the ability of these harmonics to probe the nature of DM velocity 
substructure and the extent to which it is present in the galaxy. In Sectionj^we provide rough 
estimates of the minimum number of events that would be needed to confirm the existence 
of an annual modulation for the different halo models we consider, as well as the number 
of events needed to differentiate between these models. A summary and our conclusions are 
provided in Section 


2 DM signal and its modnlation 


Since WIMPs in the galactic halo are non relativistic, u/c ~ 10“^, the amplitude for DM 
scattering with a nucleus initially at rest is usually expanded in powers of the small WIMP 
speed V and momentum transfer q = \/2nvpE^, with rriT the mass of the target nucleus and 
Er its recoil energy. The zeroth order term of the scattering amplitude in this expansion 
is u-independent and is usually the only term retained (unless it vanishes). This makes the 
angular differential WIMP-target nuclide (T) cross section dar/dcosO independent of the 
WIMP speed. However, a more useful quantity entering the analysis of experimental data is 
the differential cross section in recoil energy dar/dEji. For elastic scattering, Er = u^(l — 
cos 9)/rriT, with /iy the WIMP-nucleus reduced mass, and one finds dar/dER oc 1/u^. 

This proportionality also holds in the leading order term for inelastic scattering. 

When the scattering amplitude is independent of v, we have dor/dER = aT{E^)mT/{2^‘^v'^), 
where ariE^) is a factor with units of a cross section. Gonsidering the additional factor of 
V coming from the DM flux, the scattering rate reads 


dRx 

dER 


Ct 


o't(Er) 


pr]{v 

min 


( 2 . 1 ) 


where Ct is the target mass fraction in the detector, p is the local DM mass density, 
'ymin(ER) = JmTEji/2p^ is the minimum WIMP speed necessary to induce an elastic scat- 
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tering event with nuclear recoil energy E^, and 
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with f{v,t) the local DM velocity distribution in Earth’s rest frame. 

To interpret the outcome of direct DM detection experiments within models of particle 
DM, one typically needs to assume a specific form of the DM velocity distribution. Once this 
is specified in some reference frame R, e.g. the galactic rest frame, the velocity distribution in 
Earth’s rest frame can be obtained by the Galilean transformation f(v, t) = f-R^{v + UEs(i) + 
where /r, is the velocity distribution in R, VEsit) is Earth’s velocity with respect to 
the Sun, and usr is the Sun’s velocity in R. The time dependence of i»es( 0 is due to the 
annual rotation of Earth about the Sun. For our analysis we take v^sit) from Ref. |22j . 

In the Sun’s reference frame, DM particles that are on average at rest with respect 
to the galaxy are seen to have a preferred direction of motion that opposes r>sR- For this 
reason, the Sun experiences a constant “wind” of DM particles. The gravitational potential 
of the Sun bends the trajectories of these DM particles, acting as a gravitational lens that 
focuses the DM particles on the leeward side. As a consequence, the DM density and velocity 
distribution acquire a dependence on Earth’s relative location to the Sun (see e.g. Fig. 1 of 
|24j for a diagrammatic representation). 

The effect of GF is taken into account by replacing fR{v + VEs{t) + vsr) with /r(usr + 
VooivEsit) + '*^])) where Uoo['i^] is the velocity a WIMP had very far away from the Sun, where 
the Sun’s gravity is negligible, such that its velocity when it reaches Earth is v [231 
Ref. [23] has shown that Uoo[i’] is given by 


■n = 


-. v + ^VpcU^scr - Vocv{v • r) 

'jL + ■ r) 


(2.3) 


where Uesc = \/2 GMq7^ ss 40 km/s is the escape velocity of the Solar System at Earth’s 
position, r is the Sun-Earth distance, and f is the unit vector pointing from the Sun to Earth. 
Energy conservation ensures that 


'*'esc* 


The velocity integral in Eq. (2.2) can be written as a Fourier series. 


OO 

vivrain, t) = ao(umin) + ^ (an(t'min) cos[na;(f - to)] + 6n(t'min) sin[na;(f - to)]) , (2.4) 

n=l 


with oj = ‘I'KIjqrx and to an arbitrary phase parameter. If Earth’s orbit is assumed to be 
perfectly circular and the DM velocity distribution is isotropic, choosing to to be the time 
at which the speed of Earth with respect to the galaxy is maximum simplifies Eq. (2.4) by 
setting all = 0. Accounting for astrophysical uncertainties in the velocity of the Local 
Standard of Rest, the time at which the speed of Earth with respect to the galaxy is maximum 


occurs somewhere between May 30^“ and June 2“ |26) . 

While the coefficients of Eq. (2.4) are more easily computed, we find that it is more 
intuitive and accessible to characterize each harmonic by a single amplitude and a single 
'I’min-dependent phase: 


OO 

vivmm, t) = Ao(Uniin) + ^ A„(Umin) COs[na;(t - tn{Vmm))] , 

n=l 
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with all An ^ 0. Comparing Eq. (2.5) with Eq. (2.4) we find 


An \/®n ) 


tn = -arctan 

nuj 


an sin (nujto) — bn cos {nuto) 
an cos [ncoto) + bn sin (nuito) 
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In the next section we compare the amplitudes and phases of the first few Eourier modes of 
r]{vniin,t) for a variety of DM velocity distributions, both with and without GE. These can 
be computed analytically if one excludes the contribution of GE, otherwise the calculation 
needs to be done numerically |24| . 

When considering more than one halo component (e.g. SHM plus a stream or a DD), 
we assume that the DM consists of a single type of particle, unless otherwise noted. In 


this case the amplitude of each mode in the expansion (2.4) for the various halo compo¬ 


nents can be combined linearly, weighted by their density contribution, i.e. a„ (umin) = 
Ptot)an,iivinm) (and analogous equation for 6„(umin))) where i labels the various DM 


subcomponents and ptot = Yli Pi- Notice that the amplitude An of the expansion in Eq. (2.5) 
is not in general obtained by linearly combining the Eourier amplitudes of different DM 
subcomponents, since the phases will in general be different. If the DM is assumed to be 
composed of multiple types of particles with different masses and/or interactions, the combi¬ 
nation is not so straightforward and one must take into account the whole factor multiplying 
r/(umin) t) in Eq. (2.1). This will be relevant for the case of DM with a dissipative component 


in Section 3.3 


3 Modulation analysis for various halo models 
3.1 SHM 

In the SHM the DM velocity distribution in the galactic rest frame is described by an isotropic 
Maxwellian truncated at the galactic escape speed Uesc: 

^ 2 1^2 

/r(^) = . '1x3/2° " 1^1) ’ (3-^) 

(^^o) -Nesc 

where vq is the most probable speed in the galactic rest frame and the normalization is chosen 
so that / d^v f^{v) = 1, 


iVesc = Erf(uesc/^^o) " . (3.2) 

Vvruo 

We assume vq to be equal to the speed of the Local Standard of Rest, vq = 220 km/s [27] . 
we take the Sun’s velocity with respect to the galaxy from Ref. [22], usr = (11, 232, 7) km/s, 
and Vesc = 533 km/s [28] . 

Eig. shows the amplitudes and phases of the first (solid lines) and second (dashed 
lines) harmonic including both GE and the eccentricity of Earth’s orbit in the calculations 
(black), including the eccentricity but neglecting GE (yellow), including GE but neglecting 
the eccentricity (red), and neglecting both the eccentricity and GE (blue). Notice that, due to 
the negligible impact of the eccentricity on the first harmonic, the solid black line extensively 
overlaps the red, and the solid dark yellow line completely overlaps the blue, for both the 
amplitude and the phase. In agreement with Refs. [241125] , Eig. shows that the inclusion of 
GE causes an approximate 20 day shift in the phase of the annual harmonic for Umin ^ 100 
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Figure 1. Amplitudes (left) and phases (right) for the first (solid) and second (dashed) harmonics, 
for the SHM including the effect of GF and the eccentricity of Earth’s orbit (black), including GF but 
neglecting the eccentricity (red), including the eccentricity but neglecting GF (yellow), and neglecting 
both GF and the eccentricity (blue). The eccentricity has a negligible impact on the first harmonic, 
therefore the solid black line extensively overlaps the red, and the solid dark yellow line completely 
overlaps the blue, for both the amplitude and the phase. 


km/s and eliminates the sudden phase flips, i.e. the occurrence of a jump in the phase of a 
given harmonic by half the period. The phase flip can also be identified by the vanishing of 
the amplitude of a given harmonic. When GF is included, the amplitudes no longer vanish 
and the phases develop a softer Umin dependence. Actually, the presence of any anisotropy 
eliminates the phase flip and leads to continuous transitions in the phase. For the remainder 
of this paper we will loosely use the term phase flip to refer to both the previously defined 
jump in the phase, and the rapid, but continuous, phase transitions that may appear when 
anisotropy is present. 

The conclusion of Ref. [22] that the ratios of the amplitudes of various harmonics can 
be used to probe the anisotropy of the DM halo was based on the assumption that the DM 
velocity distribution is isotropic in the galactic frame. The existence of GF is not consistent 
with this assumption. Thus the ratios in Ref. | 22 | only hold for large Umin) where the effect of 
GF is not significant because DM particles spend little time in Sun’s gravitational potential. 
Fig. I shows the ratio of coefficients 6i and ai, defined in Eq. ( |2.4| ), with GF (solid blue 
line) and without GF (dashed purple line). One can see that without GF this ratio would 
be independent of the scattering energy and remains at a constant value of ~ 1/59, as found 
in Ref. |22j . GF significantly alters this result for 


^mm < 300 km/s. 


3.2 Sagittarius Stream 

Since the parameters governing the velocity and dispersion of the Sgr stream are not well 
known, we will follow the assumptions of [6|. In the galactic coordinate system, we take the 
mean velocity of the Sgr stream in the galactic frame to be uggr = (—65,135, —249) km/s, and 
model this component with an isotropic Maxwellian with dispersion vq = 30y^2/3 km/s ~ 25 
km/s. We model the smooth virialized component of the halo with the SHM. 
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Figure 2. Ratio |6i/ai| (see Eq. (2.4)) for the SHM, with GF (solid blue line) and without GF 
(dashed purple line). Without GF this ratio is approximately 1/59. 


Ref. [T7] recently studied the effect of the Sgr stream with a self-consistent N-body sim¬ 
ulation and found that the addition of the stream can noticeably alter //(uminj t). Specifically, 
they found four major changes. First, incorporating the Sgr stream in a realistic halo model 
with a baryonic disk produced a 10-20% increase in the direct search event rate for values 
of Umin larger than the typical relative stream speed. Additionally, Ref. [I7j found a 20-30% 
decrease in the fractional modulation amplitude defined as 

^ _ 'Hmax ~ ^min 

Vraax 4“ ^min 

where r/max('^’min) and ?/min('^’min) are the maximum and minimum of r]{vmm, t) in time. If the 
modulation is perfectly sinusoidal with a period of one year, A coincides with the amplitude 
of the sinusoid normalized by the unmodulated component of r]{vmin,t), i.e. Ai/Aq. Further¬ 
more, the phase flip of the dominant harmonic was found to occur about 10-15 km/s lower 
in Umin, and deviations of up to about 20 days were found in the phase of the modulation 
at values of Umin near the typical speed of a WIMP belonging to the stream as seen in the 
Sun’s reference frame. We begin by investigating if and how these conclusions change when 
GF is taken into account. 

The left panel of Fig. contains the relative increase in the unmodulated component 
of r]{vmin,t) when the Sgr stream is added to the SHM (Sgr-|-SHM), for stream densities 
ranging from 1% to 5% of the DM halo density. While Fig. does show an increase in the 
unmodulated component, GF appears to have no additional effect. The reason for this is 
explained in further detail in the following paragraphs. 

The right panels of Fig. [^show the fractional amplitude as defined in Eq. ( |3.3[ ) for the 
Sgr-|-SHM, for /Osgr/psHM = 0.01 (dashed) and 0.05 (solid), with (top) and without (bottom) 
GF. Here and in the following, the shaded regions between the dashed and solid lines highlight 
where lines corresponding to intermediate densities lie. As with the unmodulated component 
of the rate, the fractional modulation amplitude for Sgr-|-SHM seems to not be affected by 
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Figure 3. Left: Unmodulated component of r]{vniin,t), Aq, for Sgr+SHM, normalized by the un¬ 
modulated component for the SHM alone, for psgr/psHM = 0.01 (red), psgr/psHM = 0.03 (purple), 
and Psgr/PSHM = 0.05 (blue), with GF (solid) and without GF (dashed). Right: Fractional modula¬ 
tion amplitude, as defined in Eq. (3.3), for Sgr-|-SHM with GF (top panel) and without GF (bottom 
panel), and for stream densities ranging from psgr/psHM = 0.01 (dashed) to psgr/psHM = 0.05 (solid). 
Results are compared with the SHM alone, with GF (black) and without GF (yellow). 


GF. Small deviations from the SHM do occur at fmin ^ 200 km/s, and at values of Umin near 
the speed of the stream in Earth’s frame, regardless of GF. However, Fig. 4 of Ref. m shows 
that the uncertainty in the velocity distribution of the virialized component of the dark halo, 
specifically deviations from the assumed Maxwellian distribution, have a larger influence on 
A(umin) than the addition of the Sgr stream to the SHM. 

The effect of GF is known to increase as the relative WIMP velocity in the solar frame 
decreases. This is because slower WIMPs spend more time in the Sun’s gravitational poten¬ 
tial. The parameters of the Sgr stream chosen in this paper imply the majority of the WIMPs 
coming from the Sgr stream move at roughly u|gj. = 300 km/s in the Sun’s frame, with very 
few traveling below 250 km/s. Thus we do not find it surprising that GF has little impact on 
the unmodulated component of 'r]{vjnin, i) and the fractional modulation amplitude when the 
Sgr stream is added to the SHM. Notice that the SHM component of the Sgr-|-SHM has a 
lower average WIMP speed and a much larger dispersion with respect to the Sgr component, 
implying WIMPs from the background component of the dark halo will be more affected by 
GF than those in the stream. Figs. and show how this influences the amplitudes and 
phases of the first and second harmonics. 

Fig. 0 shows the amplitude Ai (left) and the phase ti (right) of the first harmonic. 
Without GF, Ai would experience almost no deviation from the SHM, except around the 
Umin values where the SHM has a dip due to the phase flip, and in a small region near Ugg^. 
where it decreases by at most a factor of 2. Including GF does not noticeably affect these 
results. When GF is accounted for, the flip of ti in the Sgr-|-SHM is delayed relative to the 
SHM alone, but it occurs more rapidly, causing a deviation in the phase by up to two months. 
As found in Ref. mi, when the Sgr stream is added to the smooth component of the halo 
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Figure 4. Amplitudes (left) and phases (right) of the first harmonic for the Sgr+SHM. Dashed 
lines are for a low density stream, psgr/psHM = 0.01, while solid lines are for a high density stream, 
PSgr/psHM = 0.05. Blue lines and regions include GF, while purple lines and regions neglect GF. 
Results are compared with the SHM alone, with GF (black) and without GF (yellow). 


there appears to be a significant deviation of up to two months in the phase of the annual 
modulation near Umin — I'ggj.- This effect occurs regardless of GF, and the inclusion of GF 
even appears to enhance this deviation for low density streams, resulting in an approximately 
20 day phase shift for psgr/psHM = 0.01. 

Fig. El shows the amplitude A 2 (left) and the phase t 2 (right) of the second harmonic. 
The results for the phase of the second harmonic are very similar to those of the first. Without 
GF, A 2 would deviate signihcantly from the SHM only at Umin — ^sgr around the point 
where the SHM amplitude has a dip due to a phase flip. When GF is included the dip 
vanishes but the enhancement of the amplitude at Umin — ^Sgr remains. This enhancement 
is roughly a factor of 4 for /Oggr/psHM = 0.05. Without GF, the phase would exhibit a flip 
at smaller Umin values and would experience strong deviations from the SHM of up to ~ 50 
days for 250 km/s < Umin < 350 km/s. The inclusion of GF again washes out the low Umin 
deviations, but leaves those above Umin ~ 200 km/s intact. Thus one would expect a 1% 
density stream to have a second harmonic phase similar to the SHM’s, except in the region 
Umin ~ where the stream could cause deviations of as much as 45 days. 

To summarize, including GF in Sgr+SHM calculations does not affect the unmodulated 
rate or the fractional amplitude. Furthermore, it appears to wash out characteristics that 
may allude to the potential existence of DM substructure in the phases of the first and second 
harmonics for values of Umin ^ 200 km/s. The features arising from the substructure near 
Umin ~ I’sgr’ however, are left intact. Since the effect of GF is negligible at U min > 300 km/s, 
should the Sgr stream contribute non-negligibly to the local DM density, anisotropies arising 
from the Sgr stream could be probed using the ratios of harmonics as proposed in Ref. [22] . 

3.3 Dark Disk 

We consider here two distinct types of DDs. The first could form from accretions of massive 
satellites onto the galactic disk In this scenario, the DM in the halo of the satellite 


- 9 - 












''min [km/s] v^in [km/s] 


Figure 5. Same as Fig. but for the second harmonic. 


galaxies and the DM comprising our own galaxy’s halo is expected to be of the same type, 
i.e. non-dissipative in nature. While the DM in the halo must be non-dissipative in order to 
maintain its known spatial distribution, Ref. m has shown that a subdominant portion of 
at most 5% could be dissipative. If such a component exists, it would form a DD in much 
the same way the baryonic matter dissipates energy and forms the visible disk. We refer to 
this second scenario as dissipative dark disk (DDD). Densities of DDs are expected to range 
from Pdisk/phaio = 0.2 to 1, with this ratio being strictly less than 3 and likely less than 2 

m- 

We model both the non-dissipative DD and the DDD using a truncated Maxwellian. 
Consistent with the values found in numerical simulations, we model the DD with a rotation 
velocity 50 km/s slower than the Local Standard of Rest (uiag = 50 km/s) and with velocity 
dispersion vq = 70 km/s |llj . While we present our results for a non-dissipative DD together 
with a SHM halo (DD-I-SHM), results for the DDD are obtained without a background 
halo component, as a specific particle model for all (dissipative and non-dissipative) DM 
components would be required before the velocity integrals of SHM and DDD could be 
combined. 

The left panel of Fig. shows the enhancement of the unmodulated component of 
r]{vrain,t), for a DD combined with the SHM (DD-I-SHM) relative to that of the SHM 
alone, for pdd/pshm = 0.2 (green), 0.5 (red), 1 (purple), and 2 (blue), with (solid) and with¬ 
out (dashed) GF. When GF is neglected, adding a DD to the halo can increase the unmod¬ 
ulated rate by as little as 150% for pdd/pshm = 0.2, or as much as 775% for /Odd/pshm = 2, 
but only for Umin ^ 200 km/s. GF does not appreciably change this result for the low density 
DD, but can increase the unmodulated component of the high density DD by an additional 
60%. 


The right panel of Fi g. shows the fractional modulation amplitude A(umin) for the 
DD-I-SHM, defined in Eq. (3.3). For values of Umin below 300 km/s, where one expects GF 
and a DD with small uiag to be most influential, the SHM predicts a fractional modulation 
amplitude of at most ~ 4%. Without GF, the addition of the DD to the SHM would increase 
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Figure 6. Left: Unmodulated component of ri{vynin,t), Aq, for the DD+SHM normalized by the 
unmodulated component for the SHM alone, for pdd/pshm = 0.2 (green), pdd/pshm = 0.5 (red), 
Pdd/pshm = 1 (purple), and Pdd/pshm = 2 (blue), with GF (solid) and without GF (dashed). Right: 
Fractional modulation amplitude, as defined in Eq. (3.3), for the DD+SHM with (top) and without 
(bottom) GF, and densities ranging from pdd/pshm = 0.2 (dashed) to pdd/pshm = 2 (solid). Results 
are compared with the SHM alone, with GF (black) and without GF (yellow). 


A to as much as 12%. The presence of GF reduces the influence of the DD, resulting in at 
most an 8% fractional modulation amplitude. For Umin > 250 km/s, the influence of the DD 
and GF vanish and the DD+SHM results are identical to those of the SHM. 

The left panel of Fig.j^shows the amplitude of the first harmonic. Without GF (purple), 
the addition of the DD would either increase or decrease the relative amplitude of the first 
harmonic, depending on the values of pdd/pshm and Umin- Due to GF the addition of the 
DD primarily enhances the amplitude of the first harmonic. This enhancement can be as 
large as a factor of 5 for pdo/pshm = 2. 

The right panel of Fig. shows how the DD impacts the expected phase of the first 
harmonic, with (blue) and without (purple) GF. Regardless of whether or not GF is included, 
the phase of the first harmonic of the DD+SHM looks identical to that of the SHM for 
Umin > 250 km/s. Without GF, the phase of the first harmonic would deviate from the SHM 
by nearly 6 months for Umin between 70 km/s and 200 km/s. When GF is accounted for, 
this phase difference between the DD+SHM and SHM is at most 4 months, and the range of 
Umin at which this deviation occurs is reduced to 100 km/s ^ Umin ^ 180 km/s. 

Fig. shows the effect of a DD on the second harmonic. As for the first harmonic, the 
DD affects neither the relative amplitude nor the phase for Umin > 250 km/s. Below this 
value, the relative amplitude of the second harmonic of the DD+SHM is primarily enhanced 
when GF is neglected, except in very narrow regions around Umin — 50 and Umin — 140 km/s, 
depending on poD- When GF is included, the amplitude for the DD+SHM is enhanced by 
up to a factor of 10 for pdd/pshm = 2 and a factor of 3 for pdd/pshm = 0.2, but only 
at values of Umin below 180 km/s. A slight reduction in the amplitude occurs for all DD 
densities plotted between U min values of 180 km/s and 250 km/s. The phase of the second 
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Figure 7. Amplitudes (left) and phases (right) of the first harmonic for the DD+SHM. Dashed lines 
are for a low density DD, pdd/pshm = 0.2, while solid lines are for a high density DD, pdd/pshm = 2. 
Blue lines and regions include GF, while purple lines and regions neglect GF. Results are compared 
with the SHM alone, with GF (black) and without GF (yellow). 


harmonic without GF would consistently differ from the SHM by up to 75 days for values of 
'Wmin below 250 km/s. GF slightly reduces the difference in t 2 between DD+SHM and SHM 
alone for Umin < 250 km/s. The maximum phase difference between DD+SHM and the SHM 
is roughly 50 days, but may be as little as 15 days for pdd/pshm = 0.2. 

In summary, we find that the existence of a DD with a lag speed of 50 km/s and 
dispersion of 70 km/s can significantly alter the unmodulated rate, the fractional modulation 
amplitude, and the phases of the dominant harmonics, but only for Umin ^ 250 km/s. A 
larger (smaller) uiag would increase (decrease) the Umin values at which the features associated 
with the DD appear. GF is shown to further enhance the unmodulated rate and fractional 
amplitude, but it diminishes the influence of the DD on the phases of the dominant harmonics. 
However, one should keep in mind that the relative importance of GF seen in this paper is 
dependent upon the chosen rotation velocity of the DD. 

We notice that the method of Ref. |22j for probing the nature of anisotropies in the dark 
halo with ratios of harmonics may not be useful for a DD, as the region in which deviations 
occur directly overlaps with the region made anisotropic by GF (at least in the example 
provided here). The effect of GF would be noticeably reduced should the DD rotate at a 
velocity significantly different from the rotation velocity of the Sun. This is because, in the 
Sun’s reference frame, WIMPs from the DD will be moving faster and spend less time in the 
Sun’s gravitational potential. 

We also show in Fig. |^the amplitudes (left panel) and phases (right panel) of the first 
(solid lines) and second (dashed lines) harmonics for the DDD alone, with (light blue) and 
without (purple) GF. The amplitudes for the DDD look very similar to those of the SHM, 
but they appear at much lower Umin values for the DDD lag speed we assume (50 km/s). 
This should be expected, as the DDD contains WIMPs coming from approximately the same 
direction as the dark halo but with a lower relative speed and a smaller dispersion. The most 
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Figure 8. Same as Fig. but for the second harmonic. 


notable differences occur in the phases. For all values of fmiin the phases of the first and 
second harmonic for the DDD are shifted approximately half a month earlier when compared 
to the phases of the SHM (see Fig. [^for comparison). 


4 Estimate of required number of events 


We begin by providing a rough estimate of the number of events that would be necessary 
to observe the annual modulation in the Sgr+SHM, the DD+SHM, and the SHM, using a 
very simple two-bin analysis. Let us split an annual cycle into two six-month periods, one of 
which is centered about the time of maximum of the rate and the other is centered 6 months 
later. For a fixed energy range, we denote the number of events in the two time bins 
and N-. To estimate the number of events needed to establish the existence of an annual 
modulation with a significance level corresponding to a standard deviations, we require that 

AAT = Ar+- Ar_ ^ , (4.1) 


with Wot = + N- the total number of observed events. We assume that the uncertainty 

of Nj^ and N- is \/Wot/2. Assuming the phase i is constant in the energy range considered, 
we can approximate the integrated rate as R{t) ~ iig + -Ri cos (27r(t — t)/year), where Rq is 
the unmodulated component of the rate and Ri is the modulation amplitude. For a fixed 
exposure MT, N± ~ MT[Rq/2 ± i^i/vr), where the factor of I/tt arises from integrating the 
cosine term over the temporal region defining each bin. 


Solving Eq. (4.1) for Wot in terms of Rq and Ri then yields 


f Rq 


Wot ^ ^— 1 


(4.2) 


Rq and Ri can be replaced by the integral of Aq and Ai over the energy range considered, 
since all additional constants relating Rq to Aq and Ri to Ai cancel in the ratio. Aq and 
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Figure 9. Amplitudes (left panel) and phases (right panel) of the hrst (solid) and second (dashed) 
harmonics for the DDD, with GF (light blue) and without GF (purple). 


Ai are functions of Umin and we use the relation between Umin and £'r for elastic scattering 
given below Eq. (2.1). 

Since we would ultimately like to know how distinguishable the Sgr+SHM and DD+SHM 
are from the SHM, we choose to evaluate Eq. (4.2) in an energy range where the amplitude 
and phase of the Sgr+SHM and DD+SHM deviate most strongly from the SHM. Eor the 
Sgr+SHM, this region corresponds to Umin values between 280 km/s and 300 km/s (see 
Fig. [3). This region is roughly consistent with a 1 keV bin centered at about 6.6 keV, for 
a 25 GeV DM particle scattering off xenon. Evaluating Eq. (4.2) in this region we find 
the modulation amplitude in the Sgr+SHM, for /Osgr/psHM = 0.05, requires roughly 2900a^ 
events to be detected with significance a sigma, while the amplitude in the SHM requires 
2000a^ events. The Sgr+SHM requires more events to be observed than the SHM because 
the modulation arising from the Sgr stream and the SHM are out of phase, leading to a 
reduction in the modulation amplitude as shown in Eig. [3 Eor the DD+SHM, we consider 
Umin values between 130 km/s and 150 km/s (see Eig. [ 3 . This region approximately coincides 
with a 1 keV bin centered at about 4.5 keV, for a 50 GeV DM particle scattering off xenon. 
The DD+SHM, for /Odd/pshm = 2, would require 170 q:^ events to be detected at significance 
a sigma while the SHM would require roughly 7500a^ events. The addition of the DD to 
the SHM significantly reduces the number of necessary events because the DD+SHM has a 
significantly larger modulation amplitude, as seen in Eig. [3 

An important question to ask is, if an experiment were to view the annual modulation 
with a significance of a standard deviations, would additional events be necessary in order 
to distinguish these models? It is clear from Eig. |3 that the modulation features of the 
Sgr+SHM deviate most from those in the SHM in the phase of the first harmonic. Thus the 
question to ask for the Sgr+SHM is, how many events must be observed between Umin = 280 
km/s and Umin = 300 km/s in order to distinguish a phase occurring in mid March from a 
phase occurring in late May. The first harmonic in the DD+SHM differs strongly from that 
of the SHM in both the phase and the amplitude, and thus it is important to check which 
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feature will be more easily distinguishable from the SHM. 

We will begin with a rough analysis of the number of events necessary to discern the 
difference between the amplitude of the modulation for the DD+SHM, assuming /Odd/pshm = 
2, from the SHM amplitude. Consider once again the same bin analysis previously used to 
determine the detectability of the modulation amplitude. The condition for distinguishing the 
modulation amplitude of the DD+SHM from the modulation amplitude of the SHM is simply 
that the difference between AA^dd+shm and AiVgHM niust be larger than the uncertainty, 
in the measurement of AN. This implies 

AAdd+shm — AAshm ^ oi\/Ntot ■ (4.3) 


With similar manipulations as above, one can arrive at the following condition on Wot: 

-2 


Wot ^ 


TT^a^ 


Ih\ 


4 LV-^o/dd+shm 


-^0 / shmj 


(4.4) 


Evaluating Eq. (4.4) in the energy range previously defined for the DD+SHM, we find that 


approximately 225a^ events are necessary to distinguish the amplitude of the DD+SHM 
from the SHM at a significance of a sigma. This implies that approximately 55a'^ more 
events must be detected after the annual modulation is observed in order to discriminate the 
DD+SHM from the SHM using only the amplitude of the modulation. 

We now consider how the phase of the modulation could be used to estimate the number 
of events that must be observed in order to distinguish the various models. Assume for a 
moment that an annual modulation has been detected and the number of events can be 
plotted against time to form a sinusoidal-like figure. Let us assume that one of the data points 
lies at (t*,iV), where N is the average number of events observed. Two cosine functions, one 
passing through the data point itself with phase ta, and the other passing through the upper 
end of its error bar with phase tb, can then be used to characterize the uncertainty with 
which the phase is known. The upper bound of the data point is proportional to the square 
root of the number of events in the temporal bin \/Win- We will assume events are evenly 
distributed across temporal bins, implying \/Win — \/Wot/4. Assuming AN is known, the 
conditions that by definition must be satisfied are 


AN 


cos(a;(t* — ta)) = 0 


(4.5) 


AN 


cos(w(t* - tb)) ^ 


tot • 


(4.6) 


We solve Eq. (4.5) for t* and restrict our attention to the solution that is closer in phase with 
the data. Substituting this result into Eq. (4.6) yields 

Aftot 


sin(a;Af) 


2AN 


(4.7) 


where At = ta — h- To be conservative, we choose to restrict the uncertainty in the phase to 
be at most one month, for which sin(u;Af) = sin(7r/6) = 1/2, which then implies 


Wot ^ 




4 



(4.8) 
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which coincidentally is the as same as Eq. (4.2). We note that the above analysis is only 
one sided in that it fails to account for the lower part of the error bar. The true 

uncertainty in the phase thus has a full width of two months, extending to one month to 
either side of the best-fit value. 

Since Eq. ( |4.2[ ) and Eq. ( |4.8[ ) coincide, the number of events required to distinguish 
the phases (with a two month error) of the Sgr-|-SHM or DD+SHM modulations from the 
phase of the SHM modulation are approximately the same as those required to confirm the 
existence of the modulation itself. We thus expect any experiment measuring the modulation 
in an energy range where the phases of the models significantly differ, to measure the phase 
with high enough accuracy to differentiate the SHM from the Sgr+SHM and DD-I-SHM 


5 Conclusions 

We have considered how gravitational focusing of DM due to the Sun’s gravitational potential 
would alter the time modulation of a DM signal. Previous studies have separately considered 
extracting information using a harmonic analysis [5^ and investigating how anisotropies in 
the DM halo might influence direct DM detection experiments [23H25| . The purpose of this 
paper is to unify these analyses and investigate how GE would alter the results of a harmonic 
analysis in the presence of DM velocity substructure. 

We performed our analysis on a dark halo described by the standard halo model (SHM), 
a SHM with an added DM stream as expected from the tidal disruption of the Sgr dwarf 
galaxy by the Milky Way, a SHM plus a dark disk (DD) with lag speed uiag = 50 km/s, and 
a dissipative dark disk alone (DDD) with the same lag speed. Our results for the SHM alone 
are in agreement with Ref. [25]. Additionally, the conclusion of Ref. |22) that there should 
exist ratios of the amplitudes of harmonics independent of Umin was shown to be inconsistent 
with the presence of GE at Umin ^ 300 km/s. This does not come as a surprise as the result 
of [22| assumes that the local DM halo in the galactic frame is isotropic, and GE inherently 
makes the halo anisotropic. 

Eor the Sgr stream, modeled with a velocity vsgr = (—65,135, —249) km/s in galactic 
coordinates, we found that GE is unlikely to significantly affect any DM particles coming 
from the stream, but can affect the smooth halo component, and thus can alter the relative 
contributions of the Sgr stream and the smooth halo to the velocity integral ??(umin,i)- We 
showed that by increasing the relative importance of the background halo, GE tends to 
reduce characteristic features that would otherwise be expected to appear in the phases of 
the annual and biannual harmonics from the inclusion of the Sgr stream component. In 
spite of this, GE does not eliminate the more prominent features which have the potential 
to alter the expected phase of the annual modulation by more than two months for values of 
'(^min ~ ^sgr with respect to the SHM alone, where Ugg^, is the speed of the Sgr stream in the 
Sun’s reference frame. 

Eor our DD-I-SHM analysis we considered a dark disk co-rotating with the baryonic 
disk but with a smaller rotational velocity. Since the relative velocity of DM in the DD is 
much smaller than in the SHM or Sgr stream, one would expect GE to have a much larger 
influence in this model. Indeed we showed that the inclusion of a DD has a large influence on 
the unmodulated rate, the fractional amplitude, the amplitudes of the annual and biannual 
harmonics, and the phases of the annual and biannual harmonics. However, these effects 
appear only at Umin ^ 250 km/s. 
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We also provided rough estimates of how many events should be observed in order to 
differentiate between the Sgr+SHM, DD+SHM, and SHM. We have determined that should 
an experiment measure the annual modulation in an energy range where the phase of the 
Sgr+SHM and DD+SHM differ noticeably from that of the SHM, the uncertainty in the 
measured phase will be small enough to allow for a discrimination between these models. 

Our conclusions support the idea that analyzing the harmonic series of the DM differ¬ 
ential scattering rate could potentially shed light on the distribution of DM in our galaxy. 
We have found that when DM velocity substructure is present, GF washes out some of the 
more distinctive features that would appear in the amplitudes and phases of the dominant 
harmonics were GF neglected. This is so because GF enhances the density of the low velocity 
WIMPs in the smooth halo component. However, deviations with respect to the SHM, most 
notably in the phases of the harmonics, can still persist and could provide insight into the 
astrophysical nature of DM. 
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